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Abstract 

The problems connected with Gaudin models are reviewed by analyzing model related 
to the trigonometric osp(l\2) classical ?'-matrix. The eigenvectors of the trigonometric 
osp(l\2) Gaudin Hamiltonians are found using explicitly constructed creation operators. 
The commutation relations between the creation operators and the generators of the trigono- 
metric loop superalgebra are calculated. The coordinate representation of the Bethe states is 
presented. The relation between the Bethe vectors and solutions to the Knizhnik-Zamolodchikov 
equation yields the norm of the eigenvectors. The generalized Knizhnik-Zamolodchikov 
system is discussed both in the rational and in the trigonometric case. 
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1 Introduction 



Classifying integrable systems solvable in the framework of the quantum inverse scattering 
method [[TJ, ^J, ||] by underlying dynamical symmetry algebras, one could say that the Gaudin 
models are the simplest ones being based on loop algebras and classical r-matrices. More 
sophisticated solvable models correspond to more complicated algebras: Yangians, quantum 
affme algebras, elliptic quantum groups, dynamical quantum groups, etc. 

Gaudin models |Q, ||] are related to classical r-matrices, and the density of Gaudin 
Hamiltonians 

N 

H {a) =^r ab (z a - z b ) (1.1) 

coincides with the r-matrix. Condition of their commutativity [H^ a \ H^] = is nothing else 
but the classical Yang-Baxter equation (YBE) 

[r a b(z a - z b ),r ac (z a - z c ) + r bc (z b - z c )] + [r ac (z a - z c ),r bc (z b - z c )] = , (1.2) 

where r is antisymmetric and belongs to the tensor product g ® g of a Lie algebra g, or its 
representations and the indices fix the corresponding factors in the iV-fold tensor product of this 
algebra (see Section 2). 

The Gaudin models (GM) related to classical r-matrices of simple Lie algebras were studied 
intensively (see [||, ^, [8], ^ |l(], 11, 12, [l3|, 14] and references therein). The spectrum and 



eigenfunctions were found using different methods (coordinate and algebraic Bethe Ansatz [g, g], 
separated variables [Hj, etc.). The correlation functions were evaluated for g = sl(2) by the 
Gauss factorization approach A relation to the Knizhnik-Zamolodchikov (KZ) equation of 
conformal filed theory was established [11, 12|, |13|] . 



There exists a variety of classical r-matrices with trigonometric dependence on spectral 
parameter. Although algebraic construction of integrals of motion is straightforward the 
calculation of the spectrum and the corresponding eigenstates, by means of algebraic Bethe 
Ansatz, still depends on the underlining Lie algebra and r-matrix. Moreover, many trigonometric 
r-matrices are invariant under diagonal action of the Cartan subalgebra only 

[h a + h b , r ab (z a - z b )\ = , (1.3) 

as opposed to the rational case where classical r-matrix is invariant under the action of the whole 
Lie algebra q. Hence, one can modify the Gaudin Hamiltonian ( |I~T| ) by adding a local generator 
of the Cartan subalgebra 

H a ^H^ =gh a + H^ . (1.4) 

This modification does not change the creation operators, but the Bethe equations and solutions to 
the KZ system. However, the dependence on the parameter g (a magnetic field) will be described 
by a difference equation [ i5j Wj \. 



The aim of this paper is to review problems connected with Gaudin models by analyzing 
the model related to the trigonometric osp(l\2) classical r-matrix. Results obtained here are in 
many respect similar to the ones we obtained in the case of osp(l\ 2) -invariant rational r-matrix 



[17|. However, connection of Gaudin model with magnetic field and KZ equations requires 



modification of the later by adding dynamical difference equation [ \15[ 16 ]. 

There are additional peculiarities of Gaudin models related to classical r-matrices based 
on Lie superalgebras due to .^-grading of representation spaces and operators. The study of 
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the osp(l 1 2) -invariant Gaudin model corresponding to the simplest non-trivial super-case of the 
osp(l\2) invariant r-matrix [18] started in [[T^]. The spectrum of the osp(l\2) invariant Gaudin 
Hamiltonians was given, antisymmetry property of their eigenstates was claimed, and a two 
site model was connected with some physically interesting one (a Dicke model). Let us also point 
out that recently rational and trigonometric sl(2) Gaudin models were used to describe different 
physical phenomena in metallic grains p0| ] and a condensate fragmentation of confined bosons 
[^TJ. Connection with perturbed WZNW models of conformal field theory was found in [22]. 

The creation operators used in the sl(2) Gaudin model (and similarly for sl(n) case) coin- 
cide with one of the L-matrix entry [§, g|. However, in the osp(l\2) case, as it was shown for 
rational r-matrix Jl7[], the creation operators are complicated polynomials of the two generators 
X + (X) and v + (n) of the loop superalgebra. We introduce ^-operators belonging to the Borel 
subalgebra of the trigonometric loop superalgebra £ t (osp(l\2)) by a recurrence relation. Acting 
on the lowest spin vector (bare vacuum) Bm(/^i, ■ ■ • ,jUm)^- the 5-operators generate exact 
eigenstates of the Gaudin Hamiltonians H^ a \ provided Bethe equations are imposed on parame- 
ters {fJ-j} of the states. For this reason the S-operators are sometimes refereed to as the creation 
operators and the eigenstates as the Bethe vectors, or simply i?-vectors. Furthermore, the recur- 
rence relation is solved explicitly and the commutation relations between the ^-operators and 
the generators of the loop superalgebra Ct(osp(\\2)) as well as the generators of the global su- 
peralgebra osp(l\2) C Ct(osp(\\2)) are calculated. We prove that the constructed states are 
eigenvectors of the generator of the global Cartan subalgebra h g i, but the B- vectors are not the 
lowest spin vectors anymore, as it was the case for the invariant model [p3[ O]. Analogously 



to the rational case Q17Q, a striking coincidence between the spectrum of the osp(l\2) invariant 
Gaudin Hamiltonians of spin s and the spectrum of the Hamiltonians of the sl{2) Gaudin model 
of the integer spin 2s is also confirmed in the trigonometric case. 

A connection between the instates, when the Bethe equations are not imposed on their pa- 
rameters ("off-shell Bethe states"), of the Gaudin models for simple Lie algebras to the solutions 
of the Knizhnik-Zamolodchikov equation was established in the papers [11, 12]. An explana- 
tion of this connection based on Wakimoto modules at critical level of the underlying affine 
algebra was given in [12]. An explicit form of the Bethe vectors in the coordinate represen- 



tation was given in both papers nilL 12]. The coordinate Bethe Ansatz for the instates of the 



osp(l|2) Gaudin model is obtained in our paper as well. Using commutation relations between 
the £> -operators and the transfer matrix t{X), as well as the Hamiltonians ffW, we give an 
algebraic proof of the fact that explicitly constructed instates yield a solution to the Knizhnik- 
Zamolodchikov equation corresponding to a conformal field theory. This connection permits us 
to calculate the norm of the eigenstates of the Gaudin Hamiltonians. An analogous connection 
is expected between quantum osp{\\2) spin system related to the graded Yang-Baxter equation 
[18| 24, 25|, 26] and quantum Knizhnik-Zamolodchikov equation following the lines of [27]. We 
point out possible modifications of the Gaudin Hamiltonians and corresponding modifications of 
the Knizhnik-Zamolodchikov equation, similar to the case of the sl{2) Gaudin model which was 
interpreted in [28, 29] as a quantization of the Schlesinger system for isomonodromy deforma- 
tion. 

The norm and correlation functions of the sl{2) invariant Gaudin model were evaluated in 
[f7|] using Gauss factorization of a group element and Riemann-Hilbert problem. The study of 
this problem for the trigonometric Gaudin model based on the osp(l\2) Lie superalgebra is in 
progress. However, we propose a formula for the scalar products of the Bethe states which is 
analogous to the si (2) case. 

The paper is organized as follows. In Section 2 we review main data of the quantum trigono- 
metric osp(l\2) spin system: the osp(l\2) solution to the graded Yang-Baxter equation (R- 
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matrix), monodromy matrix T(A), the transfer matrix t(X) = str T(A), its eigenvalues and 
the Bethe equations. The eigenvectors of this quantum integrable spin system can be constructed 
only by a complicated recurrence procedure [[H]] which is not given here. Nevertheless it is use- 
ful to remind the main data of the quantum integrable spin system because some characteristics 
of the corresponding Gaudin model can be obtained easily as a quasi-classical limit of these data. 
The trigonometric osp{\\2) Gaudin model and its creation operators % are discussed thor- 
oughly in Section 3. Some of the most important properties of these operators are formulated 
and demonstrated pure algebraically: antisymmetry with respect to their arguments, commuta- 
tion relations with the trigonometric loop superalgebra generators, commutation relations with 
the generating function i(A) of the Gaudin Hamiltonians, a differential identity, valid in the case 
of the Gaudin realization of the loop superalgebra. Using these properties of the B-operators we 
prove in Section 4 that acting on the lowest spin vector Q _ these operators generate eigenvectors 
of the generating function of integrals of motion, provided the Bethe equations are imposed on 
the arguments of the i?-operators. Possible modifications of the Gaudin Hamiltonians are pointed 
out, also. In particular, one of them yields Ricardson type Hamiltonian. An algebraic proof is 
given in Section 5 that constructed Bethe vectors are entering into solutions of the Knizhnik- 
Zamolodchikov equation of conformal field theory. Quasi-classical asymptotic with respect to 
a parameter of the Knizhnik-Zamolodchikov equation permits us to calculate the norm of the 
eigenstates of the Gaudin Hamiltonian. We pointed out that modification of Gaudin Hamiltoni- 
ans by adding Cartan element requires more complicated change of KZ system, as opposed to the 
rational case. Further development on possible evaluation of correlation functions is discussed 
in Conclusion. Finally, some definitions of the orthosymplectic Lie superalgebra osp(l\2) are 
given in the Appendix. 



2 Quantum osp( 1\2) superalgebra and corresponding spin system 



The quantum superalgebra U q (osp(l|2)) as a deformation of the universal enveloping algebra of 
the Lie superalgebra osp(l\2) (see Appendix), is generated by three elements h, v + , v~ [31, 32]. 
The q-deformed commutation relations between the generators are 



[h,v ± ]_ = ±v ± , [v + ,v-] + 
Its center is spanned by the q-deformed Casimir element 



q h - q- h 
q - q- 1 



(2.1) 



c 2 (q)=A(q)([h] q f+B(qy 



,1/2 



+ 



+ 2 _2 

tr , v 



+ 



+ • 



V , V 



(2.2) 

with 2A(q) = q + q~ l + \(q - q~ 1 ) 2 and B(q) = (q 1 / 2 + q~ x l 2 y 2 . The q-deformed Casimir 
element can be written also in the form 



c 2 (q) = ([/*- 1/2], 



q 1 ' 2 + q 



-1/2 V„+ 



2 _2 



v v + I q + q 



Ji-1 



-h+l\ v + v - 



(2.3) 



Like in the case of the Lie superalgbra osp(l\2) (see Appendix), there exists an element 

ci(?) = {[h - 1/2], + (q 1 / 2 + q- 1 ' 2 ) v+v-) , ( Cl (q)) 2 = c 2 (q) , (2.4) 



with a grading property [a, h] = and av^ 



-u ± ci. In the quasi-classical limit q — ► 1 the 



relations (2.1-4) are reduced to the Lie superalgebra osp(l\2) ones [33]. There is a coproduct 
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map A :U q -^U q ®U q consistent with the commutation relations (£J]) and a universal R-matrix 
1Z as an element of U q ® U q 32]. Let us write its matrix form R(q) = (p <8> p)Tl in the basis 
e\ ® ex,e\ tg> e2, ei (8) e3, . . . , e$ <8> e3 of the tensor product of two copies of the fundamental 
representation © p:U q ^ End 



R(q) 



\ 



a 

r 1 b 



1 



1 6 
1 a 

1 



\ 



(2.5) 



qJ 



here a = g — q 1 , b = q 3 / 2 — g 1 / 2 and c = (1 + g x )(g — g 1 ). Multiplying R(q) by 
the graded permutation V of V"W tg> V {X) : {V) ah - cd = {-l) p{a)pi ~ b) 5 ad 5 bc , one gets the braid 
group form of the R-matrix R(q) = R^V which has a spectral decomposition. Using the 
projectors on the irreducible representation components in the Clebsch-Gordan decomposition 
® = y( 2 ) © W 1 ) © one can represent this R-matrix in the form Q |2|] 



%) = -q~ 2 Po(q) ~ q^PM + qP2(q) 



(2.6) 



where projectors are 



Po(q) 



/0 



q l/2 + g-V2 



,3/2 



+ 9" 



-3/2 



g" 1 


g-1/2 


-1 











-1 











-1 


_ 9 l/2 





(2.7) 



A(<?) 



+ 



V 



-1 



o/ 



here f = g 1//2 + g 1//2 , and ^(g) = I — Po(q) ~ P\{q)- By the Baxterization procedure, or 
simply changing functions (A — a)/(A + a) to the trigonometric ones sinh(A — a)/ sinh(A + a) 
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in the osp(l\2)- invariant R-matrix [[T8|, 17] one gets the trigonometric ii-matrix related to the 
quantum affine algebra U q (osp(l\2)) 



- sinh(A - 2rj) sinh(A - 3g) 

smh(A + Ir}) smh(A + 6rj) 



(2.9) 



here q = e 2ri (normalized to A = R = 1). The L-operator of the quantum spin system on a 
one-dimensional lattice with N sites coincides with i?-matrix acting on a tensor product Vo <8> V a 
of auxiliary space Vq and the space of states at site a = 1,2, ... N 



Loa(X — Z a ) — i?0a(A — Z a ) 



(2.10) 



where z a is a parameter of inhomogeneity (site dependence) and Ro a (X, rf) = VR(X, rj) is the 
usual R-matrix. Corresponding monodromy matrix T is an ordered product of the L-operators 



N 



T(A; {z a }f ) = L 0N (X - z N ) . . . L 01 (X - z x ) =JJ L „(A - z a ) . (2.11) 



a=l 



The commutation relations of the T-matrix entries follow form the (Z 2 -graded) FRT-relation [Q] 

flia(A - M)ri(A)T 2 (/x) = T 2 ( fJ )T 1 (X)R 12 (X - n) . (2.12) 



Multiplying ( |2.12D by i? 12 1 and taking the super-trace over Vi (8) V2, one gets commutativity of 
the transfer matrix 



t(X) = str T(A) = ^(-l^T^A; {z a }f ) = T u - T 22 + T 33 



(2.13) 



for different values of the spectral parameter t(X)t(fi) = t(fi)t(X). 



The choice of the L-operators (2.10) corresponds to the following space of states of the 
osp(l|2)-spin system 

N m 
a=l 



The eigenvalues of the transfer matrix t(X) in this space are [18, 24, 25] 



A(A;{PiK 



n 



P°(A; R}f ) II " Mi) - 4 N } (A; R}f ) x 



i=i 



where ^ (A; {z a }f ) = Ub=i «i(A - ; j = 1, 2, 3 , 

ai(A) = sinh(A + 2t]) sinh(A + 377) , a 2 (X) = sinh(A) sinh(A + 377) , 

smh(/i — nrf) 



x Y[Si(X- Hj + r])S-i{X- Hj + 2ri) + 

M 

+ af\X;{z a }^)l[S^(X-^ + 3 V ) , 

i=i 

N 



a 3 (A) = sinh(A) sinh(A + 77) , S n (n) 



sinh^u + nrj) 



(2.14) 



(2.15) 
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Although according to ( 2.14 ) the eigenvalue has formally two sets of poles at A = fij — rj and 
A = fij — 2rj, the corresponding residues are zero due to the Bethe equations on the parameters 
{Hj} of the eigenstate [JT|, [24], |25|] 



TV 

n 



sinh(^j - Za + rj) 



M 



(2.16) 



k=l 



AJL s inh( / u j - z a -rf) 

If we take different spins l a at different sites of the lattice and the following space of states 



H = § V^ la) , 

a=l 



then the factors on the left hand side of (2.16) will be spin dependent too. 



Due to the more complicated structure of the .R-matrix (2.9) (see (2.7), (j2lj|)) than the gl(n), 
or gl(m\n) trigonometric i?-matrices, the commutation relations of the entries Tjj(A) of the T- 
matrix ( |2.1 1[ ) have more terms and construction of the eigenstates of the transfer matrix t(X) by 
the algebraic Bethe Ansatz can be done only using a complicated recurrence relation expressed 
in terms of [pOj] . It will be shown below that due to a simplification of this recurrence 

relation in the quasi-classical limit rj — > one can solve it and find the creation operators for 
the trigonometric osp{\\2) Gaudin model explicitly. Furthermore, the commutation relations 
between the creation operators and the generators of the trigonometric loop superalgebra as well 
as the generating function t(X) of the Gaudin Hamiltonians will be given explicitly, yielding the 
solution to the eigenvalue problem. 



3 osp ( 1 1 2 ) trigonometric Gaudin model 

As in the case of any simple Lie algebra, the trigonometric classical r-matrix of the orthosym- 
plectic Lie superalgebra osp(l|2) can be expressed in a pure algebraic form as an element in the 
tensor product osp{\\2) (g> osp(l\2) 



f (A) = coth(A) h®h + 



sinh(A) 



+ 



— — (e~ x v + (8) v' - e x v~ <g) v + ^j , 



sinh(A) 



(3.1) 



and it is a solution of the .^-graded classical YBE ( |1.2[ ) Q2p . This r-matrix can be decomposed 
naturally into positive and negative parts [ 34 ] 

1 



f(A) 



sinh A 



(3.2) 



sinh A V 2 



1 



~-h®h + 2X <g> X + - v ®w + + 



-A 



sinh A V 2 



1 



~-h®h + 2X + ®X +v + ®v~ 



V ' K> V 



It can also be represented in another form useful for modifications 

r(A) = coth(A) cf + 2 (X~ ® X + - X + ® X') - (v~ v + + 
here 

= h ® h + 2 (X + X" + X~ X + ) + (w + ®tr - iT ® u+) 



(3.3) 
(3.4) 
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The matrix form of f in the fundamental representation of osp(l\2) follows from (3.1) by sub- 
stituting appropriate 3x3 matrices instead of the osp{\\2) generators and taking into account 
Z 2 -graded tensor product of even and odd matrices. One can get it also as the quasi-classical 
limit rj — » from the i?-matrix {IS). Let us write explicitly the matrix form of f in the basis of 
the tensor product of two copies of the fundamental representation (g> (see Appendix) 



/ cosh (A) 



r(A) 



1 



sinh(A) 



\ 



cosh(A) 


-e" A 


2e" A 









e A 





-e" A 









2e A 


e A 


— cosh(A) 



cosh(A) / 



(3.5) 



with all the other entries of this 9x9 matrix being identically equal to zero. 

A quasi-classical limit rj —> of the FRT-relations ( 2.12| ) (i?(A;r/) = I + ryr(A) + 0(t] 2 ) 
and T(A; rj) = I + r]L(X) + 0(r] 2 )) results in a matrix form of the loop superalgebra relation, 
the so-called Sklyanin linear bracket, 



L (A), L (/i) 



r 12 (\-n), L (A)+ L (//) 



(3.6) 



Both sides of this relation have the usual commutators of even 9x9 matrices L (A) = L(A) ® ^3, 

L (m) = ^3 <8> L(/jl) and ru (A — n), where I3 is 3 x 3 unit matrix and L(X) has loop superalgebra 

2 

valued entries: 



L(A) 



h(X) -v~(\) 2X-(X) 
v+(X) v-(X) 
2X+{X) u+(A) -/j(A) 



(3.7) 



From the expression (3.2) of the classical r-matrix it is natural to assume that the L-operator has 
a triangular decomposition L± as A —> ±00 



h(g>h g i + 4X~ ®X£- 2v~ (8) u+ 
h®h g i + 4X + ®X~ + 2v + ®v 



gl 



J gi 



(3.8) 
(3.9) 



Here the first factors are generators in the fundamental representation (see Appendix) and 
the second factors are generators of a finite dimensional osp{\\2) Lie superagebra. 

The relation ([O]) is a compact matrix form of the following commutation relations between 
the generators h(X), w ± (/^), X^{v) of the trigonometric loop superalgebra £ t {osp(l\2)) 



MX) , h(n)]_ 

[h(X) , X ± (fi)]_ 



0. 



±2 



sinh(A — /i) 

_ e ( A -M) 

sinh(A — fj,) 



cosh(A - / u)X ± ( M ) - e T(A ^ ) X ± (A) 
(MA) - %)) , 
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[h(X) , v ± (fi)] 

[v+(X), v-Qjl)]. 
[«* (A), «± (/.)]. 



±1 



[X±(A), ^(/x)]. 
[X±(A), «±0*)]. 



sinh(A — jti) 
6 (A-aO 

sinh(A — fi) 

±2 
sinh(A — fi) 



(cosh(A - m)^^) " e T(A ^ ) w ± (A)) , 
(fc(A) - h(ji)) , 



sinh(A — ju) 
= [X±(A), = 0. 

(3.10) 

In order to define a dynamical system besides the algebra of observables we need to specify 
a Hamiltonian. Due to the r-matrix relation ( |3.6| ) the elements 

t(A) = i S trL 2 (A) = / i 2 (A) + 2[X+(A),X-(A)] + + [ V + (A), V -(A)]_ 

= /i 2 (A) +/i'(A) + 4X+(A)X-(A) +2u+(A)v"(A) (3.11) 

commute for different values of the spectral parameter 

t(X)t(/j) = t(/j)t(X) . (3.12) 

Thus, t(X) can be considered as a generating function of integrals of motion. The supertrace 
in ( p.l 1 ) for an even matrix {Ay} means str^4 = X^i=i( — an d we use the grading 
p(l) = p(3) = 0, p{2) = 1 (see Appendix). 

One way to show ( |3.12 ) is to notice that the commutation relation between t(X) and L(/j,) 
can be written in the form 



[t(X),L(p)] = [M(X,fi),L(ri] 



(3.13) 



where 



M(A, im) = - str f n 2 (A - /x) L (X)) - i str (r 2 2 (A - p)) , (3.14) 

and the second term is a quantum correction, which is absent if we consider the left hand side of 
(3.6), (3.1C) as Poisson brackets [35]. Also, in the sl(n) case this term does not contribute since 



is proportional to the unit matrix. Using ([33]) and pTfy it is straightforward to calculate M(A, /j,) 
explicitly 

cosh(A - n) h{p) -e x -» v~ (» 2e A ^ X~ 0) 



M(X,fi) 



sinh(A — fj) 



-(A-m) 







v (jm) 



1 

sinh 2 (A — /i) 



2e -(A-M)x+( / u) e-( A -^v+(/i) - cosh(A - n) h(jjt) 
cosh 2 (A-/x) + 1 





-2 j . S3. 15) 

cosh 2 (A-/i) + 1 

Substituting (3.15) into the equation ( 3.13| ) we obtain the commutation relations between t(X) 
and the generators of the superalgebra £ t (osp(l\2)). In particular 

4 e -(A-/i) 



4coth(A-/x) AT+(^)/i(A) 



sinh(A — n) 



X + (X)h(^+4X + (fi) 



2 e -(A-A') 

sinh(A — ju) 



( V + (A)^)- W +(/i) W +(A)) , 



(3.16) 
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[t(A), v + (fi)]_ = 2coth(X-p)v + (p)h(X)- 26 / A " ) / + (A)%)+.+(^) 

smn(A — fi) 



+ 



sinh(A — n) 



L-^-ti X + (X)v-(fi) - e- { ^- X) X+(/i)t>-(A)) . 



(3.17) 



Preserving some generality we can consider the representation space TL of the dynamical 
algebra to be a lowest spin p(X) representation of the loop superalgebra with the lowest spin 
vector S7_ 

fc(A)ft_ = p(X)Q- , u~(A)O_ = 0. (3.18) 

One can study spectrum and eigenstates of i(A) in this general representation TL. However, to 
have a physical interpretation we will use a local realization of the trigonometric superalgebra 
C t (osp(l\2)) with 

N 

n = ® v a , 

a=l 

as a tensor product of osp(l\2) representations. Then, 

JV 

h(X) = coth(A -z a )h a (3. 19) 

a=l 



N x _ z N „ A+2 

^-f smh(A - z a ) ^-f smh(A - z a ) 

a=l a=l 

where h a ,v^,X^ G End (14) are osp(l|2) generators in a representation l/ a associated with 
each site a. If in this realization one considers the limits A — > ±oo then one finds the expressions 
of the generators of the global Lie superalgebra osp(l\2) C Ct(osp(l\2)) in terms of the local 
generators 

N 

y 9 i = Y. Y ^ (3 - 22) 

a=l 

here Y = (/i.X*, w*). 

In particular, a representation of the Gaudin realization can be obtained by considering at 
each site a an irreducible representations Va of the Lie superalgebra osp(l\2) defined by a 
spin l a and a lowest spin vector uj a such that v~ uj a = and h a u a = —l a uj a . Thus, 

N " 

Q_= ® u a , and p{X) = \ (-l a ) coth(A - z a ) . (3.23) 

a=l ^ 

a=l 

It is a well-known fact in the theory of Gaudin models [||, ^] that the Gaudin Hamiltonian is 
related to the classical r-matrix ( |LT| ), (jO|) 

,(«o— Zb) r (-) _|_ p(2ft-2a) r ( + ) 



, sinh(z a -z 6 ) v ■ '"' 1 

o^ta of a 
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£ coth(z a - z b )h a h b + ^ _ - [e-^^X+X- + e^-^X-X+ 



sinh(z a — Zb) V 



and can be obtained as the residue of the operator t(X) at the point A = z a using the expansion 

t(A) =h 2 9 i + J2( Si a) , + 2 ■ , # , (3-25) 

^ Vsinh 2 (A - z a ) smh(A - z a ) J 

02(a) = h a 2 + 2 + + v+v~ — v~v+. As opposed to the rational osp(l\2)- 

invariant case [fL7|], the generating function (3.25) commutes only with one generator h g \ of the 
global superalgebra osp(l\2) 

[t(X),h gl ) = 0. (3.26) 

To construct the set of eigenstates of the generating function of integrals of motion t(X) we 
have to define appropriate creation operators. The creation operators used in the si (2) Gaudin 
model coincide with one of the L-matrix entry g]. However, in the osp(l\2) case the creation 
operators are complicated functions of the two generators of the loop superalgebra X + (X) and 

Definition 3.1 Let Bm{^\, ■ ■ ■ ,(J>m) belong to the Borel subalgebra of the osp(l\2) loop super- 
algebra Ct(osp(l\2)) such that 

.Bat Oil, . . . ,/xm) = v + (fj, 1 )B M -i(^2, ■ ■ ■ ,/iAf) 
+2X+&*) p-iy^—^W.O*, . . . , „„) ,3.27) 

w/f/i i?o = 1, B\(n) = v + (fi) and Bm = for M < 0. 77je notation B^_ 2 (fi2, ■ ■ ■ , Mm) 
means that the argument /ij is omitted. 

As we will show below, the ^-operators are such that the Bethe vectors are generated by their 
action on the lowest spin vector 0_ (3.18). To prove this result we will need some important 
properties of the i?-operators. All the properties of the creation operators Bm(v>i, • •• ,Mm) 
listed below can be demonstrated by induction method. Since the proofs are lengthy and quite 
technical we will present only two of them. 

Lemma 3.1 The creation operators Bm(h\, ■ ■ ■ ,Mm) we antisymmetric functions of their ar- 
guments 

Bm(^Xi ■ ■ ■ , MfcjMfc+lj • • • 1 Mm) = _ Bm{hi, ■ ■ ■ , Mfc+i, Mfe, • • • , hm) , 

(3.28) 

here 1 < k < M and M > 2. 
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Lemma 3.2 The commutation relations between the creation operator Bm and the generators 
v + (\), h(\), v~[X) of the loop superalgebra are given by 



M 



(X)B M = (-1)%, + (A)+2V- 



sinh(A - fij) 



(3.29) 



M 



M 



h{\)B M = B M U(A)+l>oth(A- Mi ) 



i=l 



i=l 



U +(A)B£ ) _ 1 + 2X+(A)^(-iy+ e ^ 



sinh(A — /ij) 
sinh(/ij — 



-I? 



(m) 

M-2 



(3.30) 



»-(A)B„ = (-l)"B Mt ,-(A)+^(-iy- 1 B<f 1 



+E 



e -(A-|»j) 
sinh(A — /Uj) 



(/i(A) - M/^)) 



Ai 



( sinh(A — /Ufc) sinh(/i/ c — ju,j" 
e -{X-Hi) 



Mi 
g-(w-Mj) 



sinh(^, - /i,) 



-5 



M-2 



+ «+(A)^(-l)^- 1 



sinh(A — pi) sinh(A — /Uj) 



(3.31) 



(i) (i i) 

here the upper index of BYj_ x means that the argument pj is omitted, the upper index ofB\f_ 2 

means that the parameters pi, p,j are omitted and O(J) is Heaviside function 

&(J) = 1 1 ifj> ° ■ 
U> 1 if j < 

It is useful to have explicit formulas for the commutators between the global generators and 
the i?-operators. 



Remark 3.1 The commutation relations between the generators of the global osp(l\2) ( [?. 22 ) 
and the Bm elements follow from the previous lemma 3.2. To see this we take the appropriate 



limit A — ► ±oo in (3.29), ( 3.30), ( 3.31 ). In this way we obtain 



M 



Vj t BM 



(-l) M B M v+ -2Y J {-l) 3 X + { H )B M ) _ 1 , 
i=i 

h gl B M = Bm (h g i + M) , 

M ( 
v- gl B M = {-l) M B M v~ gl + ^(-lYB^ hgl + h(pj) + 



M 



(Mj—Pk) 



smh(pj - fi k ) 



(3.32) 
(3.33) 
(3.34) 
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Lemma 3.3 The generating function of integrals of motion t(X) (3.11) has the following com- 
mutation relation with the creation operator Bm(p>i, . . . , /j>m) 

(M M 
2h(X) ^ coth (A - m) + coth ( A ~ cotn ( A ~ Mj) + M 
i=l i<j 



+ 2 E(- 1 ) J - i T^ \ v + (\)B<£ _ 1 + 2X + (\) 5 

^ v smh(A-^i) I smh(fj,i - fij) ^ 

M B (i) 

x^m(/x,)+4^(-1)- +m . y- 1 (e~( A -^)x+(A)^(/i i )-e'^- A )X+( / . t )^(A)) . 

f— ^ smn(A — ^jj V / 

(3.35) 

77ze notation used here for the operator /3m{^i) Is /3m(/Uj) = h(/j,i) + 2~^?, coth(^j — ^j). 

In the trigonometric Gaudin realization (3.17 — 19) the creation operators Bm(/J-i, ■ ■ ■ , MAf) 
have some specific analytical properties. 

Lemma 3.4 The B-operators in the Gaudin realization (3.17 — 19) satisfy an important differ- 
ential identity 

d A d ( . ., eW"^ 

7h>~ = E^T 



<9z a ^ efytj V sinh(^ - z a ) 



"a 



This identity will be a fundamental step in establishing a connection between the Bethe vec- 
tors and solutions to the KZ equation. 

The proofs of the lemmas are based on the induction method. As illustrations, we prove 
explicitly lemma 3.1 and the formula (3.29) in lemma 3.2. 

Proof of lemma 3.1. Consider M = 2 

2 e -(Mi-^2) 

B 2 (m,li2) = v + {p 1 )v + {^ 2 ) + ^-r- ( . 

smh(/xi - H2) 



Using the commutation relations ( |3.10| ) it is straightforward to check that B 2 (fJ,i, fj, 2 ) is antisym- 
metric 

^2(Ml>M2) = ~B 2 (H2, Mi) • 

Assume B^ifii, . . . ,/xjv) is antisymmetric N > 2 and for A r < M. We have to prove that 
Em^Ii ■ ■ ■ 1 Vm) is antisymmetric also. 

Consider j > 2, the antisymmetry of £?m(Mi> • • • > Mm) with respect to /Xj and /ij+i 
follows directly form the recurrence relation ( |3.27[ ) and our assumption. Namely, the terms 

B^_ 2 (n2,...,m)e- ilXl -^ ) /smh(n 1 - H ) and B^(fi 2 , mK^ 1 '"^ / smh(m - 
M(j+i)) enter with the opposite sign. 
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Therefore we only have to show the antisymmetry of Bm(hi, ■ ■ ■ , A*m) with respect to the 
interchange of //i and fi 2 . To see this we have to iterate the recurrence relation ( 3.27 ) twice and 
combine the appropriate terms 



Bm{hi, ■ ■ ■ , Hm) 



v (^2) + ^-n -X + {jii) B M - 2 (n 3 , • • • ,A*m) 



+ 2t,+ (/, 1 )X+(^)^(-l)^ 1 



sinh(^i - fj, 2 )' 



sinh(/i 2 - fij) 

J — 'J 

+ 2v+( M2 )X+0*i) Bg_ 3 ( M 3, • • • , MiVf) 

7=3 smn lR ^j) 

M -(,„-,, A M 



+ 4X+( m )X+(^)]T(-l 



1.7 e 3 ^(_!)fc+e(i~fc) 



i=3 



sinh(/ii — fij) 



k=3 



sinh(/i 2 - jttfc) 



^-4(^3, • • • ,MAf) 



(3.37) 



where i?j£_ 4 (/i3, . . . , ^m) means that the arguments //,■ and [Xk are omitted. Since v + (fi) com- 
mutes with the antisymmetry of the right hand side of (3.37) with respect to /ii and /x 2 
follows. Hence we have demonstrated the lemma. □ 
Proof of lemma 3.2. Here we prove explicitly only formula ( [3.29 ). In particular, when M = 1 
the expression ( 3.29j ) is just the anticommutator between v + (\) and v + (fi). Using the recurrence 
relations ( |3.27| ) it is straightforward to check that the formula ( |3.29 ) holds for M = 2 

v + (X)B 2 (m,fi2) = B 2 (ij>i,ij.2)v + {\) 

~ — hTT rfe-^X+CAJ-^JC+OiO) v + {^) 

smh(A — fix) V / 

+ -T7T -Je^ x -^X+(X)-e x -^X+^ 2 )) v+fa) . (3.38) 

smh(A — ^ 2 J V / 

Therefore we can proceed to demonstrate the lemma 3.3 by induction. Assume that the relation 
( j3~29l ) holds for B N , M > N > 2. Then we have to show the formula ([S~29|) is valid for M + 1. 
We use the recurrence relations (p.27|) to write 



M+l 



v+(X)B M +i = v + {\)\v + { I i 1 )B m + 2X + {ij j1 )Y j 
= -v + ^ 1 )v + (X)B M 



sinh(A — fii) 



M+l 

+ 2X+( W )E(" 1 ) 
3=2 



u + (A)S^_ 1 . 



sinh(/ii — jUj) 



e-(Mi-w) ^ 
sinh(^i — /ij) 1 



(3.39) 



Now we can substitute the expressions for v + (X)B M and v + (\)B$_ v After rearranging the 
terms in an appropriate way we have 



M+l 



v + (X)B 



M+l 



■l) M+1 B M+ iv+(\) + 2 £ . 

^ smh(A - 
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x ( e - (A -^ ) X+(A) - e-^- x) X + (^ B M ] . 

(3.40) 

This completes the proof of the lemma. □ 
The proofs of the other lemmas are analogous to the proofs we have illustrated above. They 
do not contain illuminating insights and are considerably longer than the two we have seen. Thus, 
we will omit them. 



The recurrence relation (3.27) can be solved explicitly. To be able to express the solution in 



a compact form it is useful to introduce a contraction operator d. 

M 

Definition 3.2 Let dbe a contraction operator whose action on an ordered product Y\ v + (p,j ), 

3=1 

M > 2, is given by 

M ~ l M e -(«-M*) M 



j=l k=j+l m^j,k 



(3.41) 

where cr(jk) is the parity of the permutation 

a : {l,2,...,j,j + l,...,k,...,M) (l,2,...,j,k,j + l,...,M) . 

M 

The d operator can be applied on an ordered product Yi v+ (Hj) consecutively several times, up 

3=1 

to [M/2], the integer part of M/2. 



Theorem 3.1 Explicit solution to the recurrence relation (3.2? ) is given by 



M [M/2] M M 

• • • , mm) =n v+ (^) + s m\ dm n v+ ^i) = ex P d n ■ ° ai) 

3=1 m=1 ' 3=1 3=1 

The properties of the creation operators % studied in the this Section will be fundamental 
tools in determining characteristics of the trigonometric osp(l\2) Gaudin model. Our primary 
interest is to obtain the spectrum and the eigenvectors of the generating function of integrals of 
motion t(X) ( p.lip . 



4 Spectrum and eigenstates of trigonometric osp( 1 1 2 ) Gaudin model 

With the help of the creation operators Bm it is possible to obtain the eigenvectors as well as the 
corresponding eigenvalues of the Hamiltonians of the trigonometric Gaudin model. This result 
is a direct consequence of the following theorem. 



Theorem 4.1 The lowest spin vector (3.18) is an eigenvector of the generating function of 



integrals of motion t(X) ( ]3.1 1 ) with the corresponding eigenvalue Aq(A) 



t(A)fi_ =A (A)tt_ , A (A) = p 2 (X) + p'(\) . (4.1) 
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Furthermore, the action of the B-operators ( 3.27) on the lowest spin vector J7_ yields the eigen 
vectors 

. . -,Hm) = B M (jJ>i, ■ ■ • ,A*Af) 0- , 

of f/ie i(A) operator 



w/f/i eigenvalues 
here 



A(A; {^}|£ 1 ) = ^ + 9 A y. 



y(A; {^}^i) = p(A) + ^coth(A - 



(4.2) 

(4.3) 
(4.4) 

(4.5) 



fc=i 



provided that the Bethe equations are imposed on the parameters {Hj}jL\ of the state (4.2) 



M 



/3m(/%0 = p(fij) + ^ cot M/^' - Mfe) = . 
Mi 



(4.6) 



Proof. The equation ( |4.1[ ) can be checked by a direct substitution of the definitions of the operator 
t(X) and the lowest spin vector f2_, the equations ( 3.11 ) and ( 3.18 ), respectively. 

To show the second part of the theorem, we use the equation (4.2) to express the Bethe vectors 
. . . ,p M ) 

t(\)V(ni,...,HM) =t(\) B M (fii, ■ ■ ■ , Hm) 0- • ( 4 - 7 ) 



Our next step is to use the third property of the £?-operators, the equation ( |3-3[ ), and the definition 
of the lowest spin vector the equation (3.18) in order to calculate the action of the operator 
t(X) on the Bethe vectors when the Bethe equations (46) are imposed 

t(\)B M tt- = B M t(\)n„ 



M 



+ 2p{\) ^2 coth(A - m) + 2^2 coth(A - Hi) coth(A — Hj) + M B M VL^ . 



i=l 



i<3 



(4.8) 



We can express the first term on the right hand side since we know how the operator t(A) acts on 
the vector f2_, the equation ( |4. Thus we have 



t(X) B M ^l- = A(A; {pjjfU) B M n 



(4.9) 



with 



M M 

A(A; {Hj}f=i) = A o(A) + 2p(A) ^ coth(A - Hi) + 2 ^ coth(A - /Uj) coth(A - Hj) + M , 

i=l i<j 

and we complete the proof by expressing the eigenvalue as 

M 

A(A; {pj}f =1 ) = y 2 + 8 x y , with y(X; {p 3 }f =1 ) = p(X) + £ coth(A - Mfc ) . 



k=l 



□ 
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Corollary 4.1 In the trigonometric Gaudin realization given by the equations (3.19), (3.2G), 



( 3.21 ) and ( 3.23 ) the Bethe vectors ^(p-i, ■ ■ ■ ,P-m) ( ¥■% ) are the eigenvectors of the Gaudin 
Hamiltonians (3.24) (see also /p7|/j 



# (a) *(Mi,..., 



Hm) = e m^(m^ ■ • • >MAf) 



with the eigenvalues 

E M =^2 la lb Cot H z a ~ z b) + ^2 la Coth (Mj 



N 



M 



(4.10) 



(4.11) 



6=1 



3=1 



when the Bethe equations are imposed 

M 

f3 M {Hj) = p(jij) + ^2 cotn (/"j ~ Mfc) = /X-Q coth(/ij - z a ) + } y coth(/Xj - p, k ) = . 

(4.12) 



Mi 



N 



a=l 



M 



Proof. The statement of the corollary follows from residue of the equation (4.3) at the point 
A = z a . The residue can be determined using ( 3.25| ), ( |4~4| ) and (|4~l|). □ 
Comparing the eigenvalues Ej£ ( |4.1 1[ ) of the Gaudin Hamiltonians and the Bethe equations 
( 1.12 ) with the corresponding quantities of the sl(2) Gaudin model [||, ^] we arrive to an inter- 
esting observation. 

Remark 4.1 The spectrum of the osp(l\2) trigonometric Gaudin model with the spins l a co- 
incides with the spectrum of the sl(2) trigonometric Gaudin system for the integer spins (see 
an analogous observation for partition functions of corresponding anisotropic vertex models in 



Remark 4.2 The Bethe vectors are eigenstates of the global generator h g i 



(4.13) 



As oppose to the osp(\\ 2) -invariant model [|l7|], these Bethe vectors are not the lowest spin 
vectors of the global osp(l\2) since they are not annihilated by the generator 



V gl ^(pi,...,p M ) 7^0 , 



(4.14) 



once the Bethe equations are imposed ( 4.12| ). These conclusions follow from the remark 3.1, in 
particular the equations ( 3.33 ) and ( |3.34| ), and the definition of the Bethe vectors (4.2). 

As was pointed out already in [0] for the sl{2) case, there are several modifications of the 



Hamiltonians (3.24). One of them is the Richardson's pairing-force Hamiltonian [Bq, 38, pp. 



These modifications can be formulated in the framework of the universal L-operator and r-matrix 
formalism (3^5) [f]]. 

Due to invariance of the r-matrix (|3.5[) with respect to the Cartan element 



[r(\),h® 1 + 1(8) h] = 0, heosp(l\2) 



(4.15) 
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one can add to the L-operator the element h 

L(A) -^L(X) = gh + L{X), 
preserving commutation relations ( p.6[ ). Then 

t{X) = \ str L 2 (A) = t{X) + 2g h(X) + g 2 , 



(4.16) 



(4.17) 



will have the commutativity property, i.e. i(A)i(/i) = t{^)t{X). Hence we can take t(X) to be 
the generating function of the (modified) integrals of motion 



N 



t(\) = (h gl -g) 2 + Y J 



c 2 [a 



+ 2 



~ Vsinh 2 (A - z a ) sinh(A - z a 



jfto = res t(A) =gh a + 

A = Z a 



(4.18) 
(4.19) 
(4.20) 



In this case the eigenstates are generated by the same B-operators. However, corresponding 
eigenvalues and Bethe equations are now given by 



A(A; { f i 3 }fL 1 ) = (y + g) 2 + d x y, 



here as before y(A; {/*,■} f =1 ) = Ea=i Ha) coth ( A " z <>) + Ef=i coth(A - M fc), 



?(a) 



JV M 

^ {-l a ) coth(Atj - z ) + E coth (^j - Mfc) + 5 = . 



(4.21) 

(4.22) 
(4.23) 



a=l 



The crucial step in the proof of these equations is the observation that the commutation relations 
between the operator t(A) ( |4.17 ) and the creation operators Bm are equal to the commutation 



relations (3.3) but with modified operator 0m{^]) —> 0M{fi>j)+9- To see this notice the similarity 
between the terms with v + {X)B^_ l operators and with X + (X)B ( £f} 2 operators in the lemma 
3.2 the equation ( 3.30 ) and in the lemma 3.3 the equation ( |3.3| ). 



Richardson like Hamiltonian [36, 38, pL 2C, 21] can be obtained as a coefficient in the A 



-oo expansion [17] 



t(\) = {h gl +g) 2 + 4e- 2X (^h gl -l + g)(^2e 2 ^h}j 

+ 4 *J (e e2 ^«) + Hi (e e22a ^) ) + 

Let us denote the coefficient next to the factor 4e~ 2A by H + 

H + = {h gl -l+g)\J2 e 2 ^h}j + AX+ e 2 *°X^j + 2»J ^ e 2 ^- 



(4.24) 



. (4.25) 
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This Hamiltonian is obviously not symmetric. Similar Hamiltonian can be obtained as a coeffi- 
cient in the A — > — oo expansion 

t(A) = (h gl -g) 2 + 4e 2X ^(y2e- 2 ^h}j(h gl -l-g) 

+ 4 e- 2Za *^ X* + 2 e- 2 ^ v~\ + O (e 4 *) . (4.26) 

Let us denote the coefficient next to the factor 4e 2A by H- , which is also not symmetric. Thus, we 
choose the following symmetric combination for a trigonometric generalization of the Richard- 
son Hamiltonian 

Hr = ~ (H+ + = (h gl - 1) ^cosh(2z a ) h}j +g ^sinh(2z a ) hA 

+ 2 (^(E e2 ^) + (E e ' 2 ^ + )^) 



+ 

The eigenvalues of Hr have different dependence on the quasi-momenta from the rational case 
Hr$m(vi, ■ ■ ■ ,m) = E R (M)t> M (vl, ■ ■■ ,Hm) , (4-28) 

with 

Er(M) 



(4.29) 



More complicated modifications of Gaudin models can be obtained considering quasi - clas- 
sical limit of the quantum spin system with non-periodic boundary conditions and corresponding 



M 




cosh(2/j,j) 


- ^/ a cosh(2z a )j ^ 






AI 


\ 


E sinh ( 2 ^') 


- ^/ a sinh(2z a ) g 




a=l / 



reflection equation [ 39), |40|]. The L-operator can be expressed in terms of the original one ( |3.7[ ) 

as 

L( 6GM )(A;{^}) = L(\;{z j })-L(-\;{z j }) , (4.30) 

in the case of the open chain, and it will satisfy more complicated liner brackets, defining a 
subalgebra of the loop algebra (3.6) (see also plj ] and references therein). 



Most of the trigonometric Gaudin model relations have their counterparts in the rational 
osp(l|2)-invariant case. To show this one takes a scaling limit A — ► eA, z a — > sz a , 

lim eL trig (eA; {ez a }) = L inv (X; {z a }) , (4.31) 

and in this way one reproduces known results for the osp(l\ 2) -invariant model. However, as we 
shall see in the next Section some relations of the invariant GM have quite complicated analogs in 
the trigonometric case (a generalization of KZ system to include a "magnetic field" parameter g, 
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requires a difference dynamical equation [15, 16]). Also the modified L-operator ( 4.16 ) requires 
to scale the parameter g — > g/e. 

Another modification can be obtained by performing the similarity transformation on the r- 
matrix d3.l| ) by the tensor square of the element exp(iX + ). Then the sacling limit A — ► eX, 
t — > £,/2e results in a modified r-matrix 



f (A) = -2- + £ (h ® - X+ ® h - v + ® u+) . 
A 



(4.32) 



The loop superalgebra will be modified, as well as corresponding Hamiltonians (1.1). Similarly, 
the algebraic Bethe Ansatz will require changes, although the Bethe equations and the spectrum 
will be the same as in the osp{\\ 2) -invariant case (see the sl(2) case in [42]). 

The expression of the eigenvectors of a solvable model in terms of local variables param- 
eterized by sites of the chain or by space coordinates, is known as coordinate Bethe Ansatz 
[g]. The coordinate representation of the Bethe vectors gives explicitly analytical dependence on 
the parameters {^i}f ! and {z a }\ useful in a relation to the Knizhnik-Zamolodchikov equation 
(Section 5). Using the Gaudin realization ( 3.19 ), ( 3.2CQ , ( 3.21 ) of the generators 



N 



-(/*) = £ 



{_L Za 



N 



sinh(/i - Za) 



v+, X+(v) = J2 



oft Za 



sinh(^ — z a 



and the definition of the creation operators ( p.42D , one can get the coordinate representation of 
the B-operators: 



N 



)«\[<p{{^}f a \za) , (4-33) 



a=l 



where the first sum is taken over ordered partitions it of the set (1,2,..., M) into subsets fZ a , 
a = 1, 2, . . . , N, including empty subsets with the constraints 

\JlC a = (1,2,..., M) , /C o p|£ 6 = for b. 



The corresponding subset of quasimomenta 



(a) _ (a) _ 



(a) 

Ka\ 



^i\K.„ \ 'i3m € /C a J , 



where | K a | is the cardinality of the subset JC a , and jk < jk+i, entering into the coordinate wave 
function 



i p({u m }P;z)= £ (-1) 



p{a) _ 



AW" 



7(2) 



A( 2 )~ 



7(3) 



sinh(i/ CT(1) - u ai 2)) sinh(i/ CT(2 ) - i/ ff(3) ) sinh(i/ CT( | X: |) 



Due to the alternative sum over permutations a € «5|^| this function is antisymmetric with respect 
to the quasi-momenta. Finally the first factor in ( 4.33| 



" " " V a M 



One 



means that for j rn G K a corresponding indices of u+ are equal to a so that u+ = 

can collect these operators into product n^Li i v t ) > consequently we have an extra sign factor 

(_1)pW_ 
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This coordinate representation is similar to the representations obtained in [O, [l2|, |l3|] for the 
Gaudin models related to the simple Lie algebras (see also [fl3|]). The .^-grading of superalgebra 
results in extra signs, while the complicated structure of the Bm -operators (for the si (2) Gaudin 
model they are just products of i^-operators B\(nj) = X + {nj)) is connected with the fact that 
(u+) 2 = Xj~ 7^ while for j / k and v£ anticommute. 



5 Solutions to the Knizhnik-Zamolodchikov equation 

Correlation functions ip(zi, . . . , z^) of a two dimensional conformal field theory satisfy the 
Knizhnik-Zamolodchikov equation p4|] 



(5.1) 



Kd^lpiz!,...^]*) = H^lfj(z 1 ,. • • ,Zn) 



where (a = 1, . . . , N) are the Gaudin Hamiltonians ( |3 . 24 ) and ip(z\ , . . . , zn) is a function 



N 



r(la) 



of N complex variables with its values in a tensor product H = ® Va 

a=l 

A relation between the Bethe vectors of the Gaudin model related to simple Lie algebras 



and the solutions to the Knizhnik-Zamolodchikov equation is well known for sometime [|li|,|120. 
Approach used here to obtain solutions to the Knizhnik-Zamolodchikov equation corresponding 
to conformal field theory and Lie superalgebra osp(l\2) starting from B-vectors ( |4.2| ) is based on 



[11] 



A solution in question is represented as a contour integral over the variables fj,\, . 



i Mat 



V>Oi, ...z N ) 



(j>(j2\z)V(j2\z) dm... dfi M , 



(5.2) 



where an integrating factor 4>(j2\z) is a scalar function 

M N 



M N 



4>{p\z) = fjsinh(/xj - JJsinh( 



Za ~ Z b ) 



n n smh (/ u 



k Z c 



(5.3) 



i<3 



a<b 



\k=lc=l 



and ^ (p\ z) is a Bethe vector (^) where the corresponding Bethe equations are not imposed. 

As a first step in the proof that ip(zi, . . . zn) given by (5^) is a solution of (5T) we differen- 
tiate the product <p^f with respect to z a and obtain 



d Za (0*) = d Za (0) * + <f>d Za (*) . 
Using ( ^3| ) the first term on the right hand side can be calculated explicitly 



(5.4) 



k3. 



N 



M 



^ la k COth(z a ~ Zb) - ^ l a COth(z a - flj) 



6=1 



3=1 



E 



(a), 
M 1 



(5.5) 



/ 



Furthermore, taking a residue of ( |3.3[ ) at A = z a we have 

M -(*»-^) 



3=1 



sinh(2: a - jij 



(5.6) 
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where 



A + 2X+^(-l) fc + ^ 



e ^ 



Mi 



-5 



M-2 



(5.7) 



Hence (5.4) can be written as 



K d 2a (^)=ffW (^) + 0^(-l)^ 



Moreover, from (5.3) we also have 



/ 



k<9, 



iV 



smh(jUj - z a 



-^(A*i)* W,a) + «^*.(*) • (5-8) 



\ 



> ,(-^a) COth(/ij - Z a ) + COth(/ij - /ifc) 



a=l 



fe=l 



<j) = P M ([J>j)<t> , 



(5-9) 



and from the lemma 3.4 follows 



M 



sinh(/x,- - 2T ) 



/ 



(5.10) 



Thus, using ( p.9| ) and ( p.lO| ), we can combine the last two terms in ( |5.8| ) into a sum of first oder 
derivatives in fij 



M 



Kd Za (^) = #w (0*) + 



sinh( / u i - 



j . 



(5-11) 



A closed contour integration of <fity with respect to fj,\, . . . , [am will cancel the contribution 
from the terms under the sum in (|5.11| ) and therefore ip(zi, . . . zn) given by (5.2) satisfies the 
Knizhnik-Zamolodchikov equation. 

Conjugated Bethe vectors (Bm^-)* are entering into the solution ip(zi, ■ ■ ■ zn) of the dual 
Knizhnik-Zamolodchikov equation 



- n—i>{zi, . . . ,z N ) = ip(z 1 , . . . ,Z N ) . 



(5.12) 



The scalar product yp(zi, . . . , z^) , i/j(zi, . . . , zn)j does not depend on {zj}^ and its quasi- 
classical limit k — > gives the norm of the Bethe vectors due to the fact that the stationary points 
of the contour integrals for k —> are solutions to the Bethe equations [13] 

N M 
y^(-^q) COth(/ij - Z a ) + ^ COth(/ij - /Xfc) = , (5.13) 



dS 



3 a=l 



k=l 
3+k 



N 



K In cf> = l a l b In (sinh(z a - z b )) 



a<b 



M N M 

+ ^ ln (sinh(/ij - /ij-)) - y~] ^ / a In (sinh(z a - fij)) 
i<j a=l j=l 



(5-14) 
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According to the remark in the end of Section 4 analytical properties of the Bethe vectors of the 
trigonometric osp{\\2) Gaudin model coincide with the analytical properties of the trigonometric 



sl(2) Gaudin model. Thus, the expression for the norm of the Bethe vectors ^ (4.2) obtained as 
the first term in the asymptotic expansion k —> coincides also 



d 2 s 



= det 



la 



d 2 s 



N 

Y- 



M 



(5.15) 



E 



1 



d 2 s 



sinh (/ij 



djij sinh 2 (^- — /Ufc) ' 



(5.16) 



for j / k. 

Finally we notice that the modification of the Gaudin Hamiltonians we discussed at the end of 
the previous Section, can be easily transfered to the corresponding modification of the Knizhnik- 
Zamolodchikov equations. The modification ( 4.16| ) for the s/(2)-invariant Gaudin model was 
studied in [29] as a quantization of the Schlesinger system (see also [28]). This modification is 



related with extra factor in the integrating scalar function (5.3) 



S 3 

exp I — 

K 



j = o, 1 , 



(5.17) 



where S = 5(5.14) and 



(5.18) 



M N 

Si = S + g ^ m - g ^ l a z a , 

j=\ a=l 

correspond to the modification ( 4.16| ). 

Moreover, following the lines of $Sft, one can try to extend the connection between the KZ 



equation and the Guadin model based on the modified L-operator ( |4.16| ) by extending the KZ 
system to include an equation of the form 



9 n 

^7 n Rich 

dg 



ip = 0. 



(5.19) 



However, such a sttaightforward generalization has failed in the trigonometric case (see below). 
We can comment on the extension in the rational case [|l7|] as a scaling limit of the trigonometric 



Gaudin model ( |4.3ip . The equations of original KZ system are defined by mutually commuting 

cf(a,b) 



differential operators (see (3.4)) 





dz a 



H 



(«) 



d 



dz a 



gh a - Yj 



Zb 



(5.20) 



The operator 



d_ 

dg 



Hr 



d 



N 



id i 



Yzghg 



a=l 



2.9 



{c 2 {gl) - hgi(h g i ~ 1)) 



(5.21) 



is commuting pairwise with the operators V a . Thus in the rational case the KZ system can be 
generalized to include the operator V 9 . 
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To prove that the solution to the modified KZ system with <$>\ given by ( 5.17 ) and ( 5.18|) is 
a solution to the generalized KZ system we have to extend the trigonometric KZ equations (5.1) 
with modified Hamiltonians gh a + along the lines of [15, 16]. A difference equation must 
be introduced 



K(zi, ... ,z N ;g) ip(zi,. . . z N ;g) = ip(zi,. ..z N ;g- 2k) 



instead of (5.21). The operator K is defined on the space Tt 
K(zi, ...,z N ;g) 



exp ^-2^z a /i a ^ P(g;h g i,v^,v gl ) 



(5.22) 



(5.23) 



where the operator P depends on the global generators of the subalgebra osp(l\2) C £t, and 
is constructed form the extremal projector p(h, v + ,v~) by a shift of the Cartan generator (see 
Appendix). We introduce only one K operator since the rank of osp(l\2) Lie superalgebra is 
one. In general case, of simple Lie superalgebra of rank r, one has to consider a set of Kk, 
k = l,...,r(seeft&M 



6 Conclusion 

By analyzing the model related to the trigonometric osp(l\2) classical r-matrix the algebraic 
Bethe Ansatz approach to the Gaudin models is reviewed. The results presented in this paper 



are is some sense analogous to the ones we obtained for the osp(l| 2) -invariant model [17]. In 
particular, a striking similarity between some of the most fundamental characteristics of this sys- 
tem and the si (2) trigonometric Gaudin model was confirmed. Although explicitly constructed 
creation operators E>m ( |3.42[ ) of the Bethe vectors are complicated polynomials of the L-operator 
entries v + (X) and X + (X), the coordinate form of the eigenfunctions differs only in signs from 
the corresponding states in the case of si (2) trigonometric model, being antisymmetric functions 
of the quais-momenta. Moreover, the eigenvalues and the Bethe equations coincide, provided 
that the sl(2) Gaudin model with integer spins is considered. Analogously, the KZ equations 
based on both trigonometric models and for the nontrivial magnetic field g require extension of 
the system of equations by the dynamical difference equation. 

Let us point out that by the method presented in this paper one can construct explicitly cre- 
ation operators of the Gaudin models related to trigonometric Izergin-Korepin r-matrix [[| |30| ] 

(2) 

corresponding to the twisted affine algebra A 2 . Similarly to the simple Lie algebra case solu- 
tions to the Knizhnik-Zamolodchikov equation were constructed from the Bethe vectors using 
algebraic properties of the creation operators Bm and the Gaudin realization of the loop superal- 
gebra Ct{osp{l\2)) . This interplay between the Gaudin model and the Knizhnik-Zamolodchikov 
equation enabled us to determine the norm of eigenfunctions of the Gaudin Hamiltonians 



V(tii,...li M ;{za}T) ||=det 



The difficult problem of correlation function calculation for general Bethe vectors 

c ({vj}?\ {Mf ; {Afelf ) = (n- , b* m {u 1} . . . u M ) JJ h(\ k ) b m (jh, . . . mm)o_^ 

was solved nicely for the s/(2)-invariant Gaudin model in [0] using the Gauss factorization of 
the loop algebra group element and the appropriate Riemann-Hilbert problem. Although the 
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corresponding factorization is known even for the quantum superalgebra U q (osp(l\2)) [ fi3| ] the 
final expression of the correlation functions is difficult to obtain due to the complicated structure 
of the creation operators Bm(hi, ■ ■ • Mm) = Poly(w + , X + ) ( p.42 ). The study of this problem is 
in progress and the following expression for the scalar product of the Bethe states is conjectured 
(see [§) 

, B* M {y 1 ,...v M )B M {n u ...n M )toJ) = Y, (-l) p(a) det>[ CT , 

<t£Sm 

where the sum is over symmetric group Saj and M x M matrix M° is given by 

■Ml* = -r—n r (p(fj>j) - p{v a (j))) - > -r— 77 \ ■ r? \ ' 

JJ smh^j-i/^)) ^ smh(/Uj - fj, k ) smh(u aij) - u aik) ) 



smh(nj - fi k ) sinh(z^ (j) - v a(k) ) 



M° k = — r— ^— — 7 , lor././,- 1.2 \I. 
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A Appendix: Orthosymplectic Lie superalgebra osp(l|2) 



The rank of the orthosymplectic Lie algebra osp(l\2) is one and its dimension is five [33]. The 
three even generators are h, X + , X~ and the two odd generators are v + , v~. The (graded) com- 
mutation relations of the generators are 



[h,X ± ]_ = ±2X ± , [X+,X~]_=h, 

[h,v ± ]_=±v ± , [v+,v-] + = -h, 
[XT, v ±]_ = V T t [ u ± j w ±] + = ±2X± , 



(A.l) 



The Casimir element is 

c 2 = h 2 + 2 (X + X~ + X~X+) + (v + v~ - v-v + ) 

= h 2 -h + 4.X+X- + 2v + v' . (A.l) 

It is interesting to point out the existance of a "square root" of this element 

Cl = h + 2v + v~ --, ( Cl ) 2 = C2 + i, (A.3) 

with a grading property [c±, X^] = 0, [ci, h] = and civ^ = —v c\. The finite dimensional 
irreducible representations V^' of the osp(l\2) Lie superalgebra are parameterized by an integer 
I, so that their dimensions 21 + 1 and the values of the Casimir element ( [\.2[ ) C2 = 1(1 + 1) 
coincide with the same characteristics of the integer spin I irreducible representations of sl(2). 
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The fundamental irreducible representation of osp(l\2) is three dimensional. We choose 
a grading of the basis vectors ej; j = 1, 2, 3 to be (0, 1, 0). Explicitly we have 



1 




















-1 











-1 











1 






(«+)* 



together with X ± = ±(u ± ) 2 . The matrix v + in the representation V® has 21 non-zero elements 
on the second upper diagonal only, and these elements are 

{(^ + )ii+i} = (A VT, VT^T, y/2, . . . , V% y/l) , 3 = 1, 2, . . . , 21. (A.4) 



The extermal projector [ |46| ] for osp(l\2) (on the lowes weight vectors): 

= p s (/ l , V + ,^)po(/i,X+,X-) , (A.5) 

here po(h, X + , X~ ) is the usual sl(2) extermal projector. There is no such factorization property 
for the extermal projector of the quantum superalgebra U q (osp{l\2)) pT[|. 
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